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ABSTRACT 


Low energy impact may be potentially dangerous for many 
highly optimized stiff" structures. Impact by foreign objects such 
as birds, ice and runways stones or dropping of tools occur frequently 
and the resulting damage and stress concentrations may be unacceptable 
from a designer's standpoint. The present work Is concerned with the 
barely visible, yet potentially dangerous dents due to impact of foreign 
objects on the Advanced Launch System (ALS) structure. Of particular 
interest is the computation of the maximum peak impact force for a 
given impactor mass and initial velocity. The theoretical impact forces 
will be compared with the experimental dropweight results for the ALS 
face sheets alone as well as the ALS honeycomb sandwich panels. 

1. INTRODUCTION. 

One of the earliest work on the low velocity impact on composite 
sandwich structures was performed by Rhodes (1974,. 1978) and:Rhodes et al (1979) 
They showed that the honeycomb core considerably reduced the area 
of delamination damage as compared with the unsupported graphite-epoxy 
laminates, and the cripping of the core allows high bending stresses in the 
face sheets. Oplinger and Slepetz (1975) found from the dropweight 
experiments that the graphite sandwich panels exhibit marked damage 
due to nominal impact energy levels as low as 2 ft-lb. due to low 
strain to fracture of graphite material and low compressive crushing 
strength of the honeycomb core. They analysed the sandwich panel 
as a plate resting on an elastic foundation. Sharma (1981) measured the 
preload and impact energy combination necessary to cause catastrophic 
failure of graphite/epoxy sandwich structures and examined the residual 
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strength of the specimens. Further, 't Hart (1981) examined the effect 
of impact damage on the tension-compression fatigue properties of 
Carbon/Epoxy Sandwich Panels and found that significant damage may occur 
at low impact energy. Labor and Bhatia (1980) examined the impact resistance 
of hybrid a«d graphite layupsand investigated the effects of core densities 
of the sandwich structure. Gottesman et al . (1987) and Bass (1986) 
presented experimental and analytical results on the strength of sandwich 
struoturesdue to low velocity impact. Further studies on the impact 
resistance and damage tolerance of sandwich plates were reported by 
Bernard (1987) and Bernard and Lagace (1987). They found that damage in the 
impacted facesheets was primarily delaminations with the largest delamination 
occurring between the bottom two plies (5*^ and 6^plies) in the top 
facesheet; and debonding of the top facesheet from the adhesive layers 
was more pronounced in stiffer core. Recent work on the instrumented 
impact testing of composite sandwich panels was performed by Shih and 
Jang (1979). They found that the impact resistance was mainly controlled 
by the facesheets and relatively independent of the density of the 
poly-vinyl -chloride (PVC) foam core, provided the facesheet material is 
tough enough; however, for less tough facesheets, the impact failure becomes 
foam core dominated rather than facesheet dominated. In particular, the 
macroscopic and microscopic failure modes and energy absorbing characteristics 
of these sandwich panels were examined by Shih and Jang (1979). 
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The NASA Langley researchers have done extensive experimental and 
analytical studies on the low energy impact resistance of graphite/epoxy 
plates. Bostaph and Elber (1982) presented. : static indentation tests 
on graphite epoxy plates and the results are compared with the theoretically 
predicted plate stiffness and maximum strain energy at failure. Bostaph 
(1984) from the US Army showed that toughened materials were able to 
delay insipient delamination, but not significantly, in quasi-static 
indentation tests. Many other experiments on impact resistance of 
composite plates have been performed by numerous investigators (see 
for example Ciarns and Lagace 1988, Sjoblom et al. 1988). For brievity, 
the remaining part of this introduction will focus on the"theoretical" 
studies on impact of composite plates. 

One of the earliest theoretical investigations on the impact of 
isotropic-homogeneous plates was presented by Eringen (1953) and 
Timoshenko in 1913 as reported in Timoshenko and Goodier (1970). 

An excellent survey of the historical contributions of various authors 
was presented by Greszczuk in 1982. Based on the Hertzian contact 
parameters for transversely isotropic laminates (Conway 1956), Greszczuk 
presented theoretical and experimental "peak" maximum impact force for 
quasi-siotropic circular plates under a spherical impactor at the plate 
center. Further investigations on the theoretical prediction of low 
velocity impact force and the delamination growth analysis in quasi-isotropic 
laminates were reported by Shivakumar and Elber (1984) and Shivakumar 
et al. (1985a, b). The present work follows closely the method presented 
by Shivakumar (1985b) which includes the finite deflection effects. 
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2. Impact Analysis of Graphite-Epoxy Plates 

The energy-balance mthod is used to obtain the peak impact force 
for circular or square graphite-epoxy plates under a concentrated 
at the plate center. This method was used by Greszczuk(1975, 1981,1982), 
Greszczuk and Chao (1977) and Shivakumar et al. (1985b) and the theoretical 
peak impact force was computed based on the equation. 


0/2)(Mj)(Vj) 2 = (l/2)(K b )W 2 + (l/4)l<y + (2/5) P 5/3 /n 2/3 

In the above expression, Mj and .Vj are the mass and velocity of the 
impactor, and are the bending and membrane stiffness of the 

plate, W fs the deflection at the center of the plate and the 
impact force P is proportional to the ralatlve displacement* raised 
3/2 power according to Hertz law (Conway 1956, Willis 1966 and 
Timoshenko and Goodier 1970) such that, 

P- 

The various coefficients of Hertz law are computed for the present 
graphite-epoxy plate and they can be found in Appendix B. From 
finite deflection plate theory, the impact force can also be written 
as (Shivakumar et al . 1985b, Vol'mlr 1967 and Timoshenko and 
Woinowsky-Krieger 1959), 


P=Kb«* 


Although the above formulae are Intended for isotropic-homogeneous 
plates, they can also be applied to transversely Isotropic material 
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(see Dahan and Zarka 1977) and quasi -isotropic materials (Agarwal and 
Broutman 1980 and Greszczuk 1982), provided one can obtain the average 
Young's modulus and average Poisson's ratio as described in Appendix A. 
The energy-balance can be written in a non-dimensional form by dividing 

q 

through by E.h (the total thickness of the composite plate is h), 

a V 


e 0 = (1/2) k. (W/h) 2 + (1/4) k (W/h) 4 + (2/5) p 5/3 /(n*) 2/3 

(4) 


where, 


e„ • 0/2)M,V, 2 /(E av h 3 ) 


k b = V< E av h > 

k m V' E av 

P ' P/(E av h2) = k b (W/h) + k m ( “ /h)3 
n* ■ n/(E av h ,/2 ) 


(5) 


The value of n for the contact between a steel punch and a graphite-epoxy 
plate is given by eqn. BIO in Appendix B, using the average value of 
the composite plate, 

E av = 7.0 xlO 6 psi, V av = 0.30 , n* = 0.567400142 

(6) 

so that the energy balance equation becomes, 

e o = (1/2) k b (w/h)2 + <1/4 ^ k m (W/h)4 + °- 583624746 P 5/3 

(7) 
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From Shivakumar et al. 1985b, there are four possible boundary conditions 

for the axisymmetric deflection of a circular plate under a concentrated 

load at the plate center. They are, using h = 0.081 In. , a * 1.50 in. 

based on tee grverning equilibrium and compatibiltiy eqns described In Appendix D. 

(i) Clamped, In-Plane Immovable, p = (616.588 lb)(W/h) + (274.273 lb)(W/h) 3 
k b = V^ E av h ^ = 4 - 603066 (h/a) 2 = 0.01342254 

k m =K m h/E av = 2.0479915 (h/a) 2 = 0.005971943 

(ii) Clamped, In-Plane Movable, P = (616.588 1 b) (W/h) + (124.012 1 b) (W/h) 3 

k b = V (E av h) = °- 01 342254 

(8a,b,c,d) 

k m = = °- 9259941 3 (h/a) 2 = 0.002700199 

(Hi) Simply Supported, In-plane Immovable, 

P = (242.846 lb)(W/h) + (347.017 1b)(W/h)3 

k b = K b^ E av h ^ = ^ 81 3329093 (h/a) 2 = 0.005287668 
k m = y^av = 2.59117215 (h/a) 2 = 0.007555858 

(iv) Simply Supported, In-plane Movable, 

P « (242.846 1 b) (W/h) + (66.09859 lb) (W/h) 3 
k fa = 0.005287668 

k m = y/ E av = 0-493558614 (h/a) 2 = 0.001439217 
The above k b and k m values for all four types of boundary conditions 
for a circular plate agree with those reported by Vol'mir (1967 section 43). 
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The bending stiffness « b should be replaced by K^c s using thick plate 
theory which takes into account transverse shear effect (Lukasiewicz 1976, 
1979 and Shivakumar et al . 1985b) where. 


V K s * MV’lH 1 - 4 W E av>H (4/3 > +1 "< 4 / a contact>] 


Using 

J = 0.28, 

rz * 

G,„ = 0.5959 xlO 6 psi 
zr 

E flV = 7.0x1 0 6 psi 

ln(a/a contact ) = ln(2a/h) = 3.6119184 


( 10 ) 


one obtains, 

K b /K $ = 1.0680294 k b (E av /G Zf ) = 12.54607 k fa (11) 

Thus, we have, for each of the four boundary conditions, 

(i) Clamped Immovable, or Movable, c = 0.855871 . k.c = 0.01148796 

(ii) s (12) 

(iii) or (iv) Simply Supported, Iirmovable or Movable, c $ = 0.9377876 

k b c $ = 0.0049587 
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For a simply supported square plate under a concentrated load P at 
the plate center, the load deflection relation is (see Timoshenko and 
Woinowsky-Krieger 1959, section 34 and Ugural 1981 section 3.3), 


W = 0.01160 Pb 2 /D 


(14) 


where b is the length of the side of the square plate and D is the 

flexural rigidity and W is the deflection at the plate center. This 
equation can be re-arranged as, letting ^ av = 0.30. 


P = 86.20689(hD/b 2 ) (W/h) = 7.894403(E a h 4 /*> 2 ) (w/h) 

a V 

or (letting b=diameter of the circular plate = 3.00 in.), 
k b = K b /(E av h) = 7.894403 (h/b) 2 = 0.00575502 

P = (264.3107 pounds) (W/h) 


05) 


(16) 


Finally, for a clamped square plate under a concentrated load at the 
center of the plate, one obtains (see Ugural 1981, section 3.12), 

W = 0.005592 Pb 2 /D > let b=2a = 3.00 in., 

P = 184.327 (hD/b 2 ) (W/h) = 16.87984 (E a h 4 /b 2 ) (W/h) 

3 V 

k b = y (E av h) = 16.87984 (h/b) 2 = 0.012305 * 

p = (565.150 pounds) (W/h) 
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Comparing the impact force on a simply supported circular plate with 
radius a and a simply supported square plate with side b=2a such that the 
circular plate is just inscribed inside the square plate, the impact 
force on the circular plate is smaller than that of the square plate. 

That is, the deflection of the circular plate is larger than that of the 
square plate for the same impact concentrated force. 

P = (242.846 lb)(W/h) + ... s.s.- circular plate, radius a 

, 08 ) 

P = (264.310 lb) (W/h) + ... s .s. square plate , side =2a 

(264.31/242.846= 1.08838) 

This somewhat unusual result arises from the fact that the deflection 
mode shape for a square plate Is quite different from that of the circular 
plate. Further, the reaction forces at the four corners of the plate tend 
to produce a convax upward deflections under the applied downward 
concentrated force P. On the other hand, for clamped plates, 

P - (616.588 lb) (W/h) + ... clamped circular plate, radius a 

, (19) 

P ~ (565.150 lb)(W/h) +... clamped square plate, side =2a 

(61 6.588/565. 1 5=1 .091 01 ) 

the effects of the reaction forces at the four corners are much less 
apparent due to the reaction moments along the four edges of the square plate. 
Since the area of the circular plate is *a 2 and the area of the square 
plate with side =2a is 4a 2 , the ratio of the area is 1.27323. The impact 
force of the circular plate is larger than that of the square plate. 

Since the experimental data lies somewhere between simply supported and 
clamped conditions, it is expacted that the impact forces for the circular 
Plates are approximately the same as the square plate. 
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3. Impact Analysis of ALS Honeycomb Sandwich Plates 


The ALS honeycomb sandwich structure is composed of a top and bottom 
graphite/ epoxy facesheets and a honeycomb core. The core is made of 
HFT-3/16-2.0 glass-phenolic material and the core material was manufactured 
by Hexcel Inc. The facesheets are identical to those examined in 
section two and each sheet is composed of 16-layer quasi -isotropic 
laminate. The material property of the honeycomb core is, 


E 


core 


= 1 7000 psi , 


G xz (L-straight direction) = 15000 psi 



(W-non-straight direction) = 5000 psi 


( 20 ) 


The above core material properties are higher-tban those for the Nomex 
honeycomb core reported by Bernard and Lagace (1987). The sandwich 
structure is being analysed as a plate resting on an elastic foundation. 

The flexural rigidity of the top sheet is, 

D = E av h3/ L120 -^v 2 )1 = 340.6673 lb-in (21) 

The elastic foundation is assumed to be linear elastic so that the 
force is proportional to the displacement and the core stiffness is, 

k = E /h = 17-00 psi/ 1.370 in = 12408 lb/in 3 (22) 

core core 

where the core thickness is 1.370 and the total thickness is 
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h total * h core + 2h - ’- 370 + (2X0-081 In) = 1.532 in (23) 


The characteristic length is 


J? = (D/k) 1/4 = 0.40705 in ( 24 ) 

Thus, from Timoshenko and Woinowsky-Krieger (1959) and Hui (1986), 
the concentrated force is related to the deflection (directly at the 
location of the application of the force) 

P = (8Dh/.* 2 ) (W/h) = (1332.324 lb) (W/h) (25) 


in the above, the top facesheet is assumed to be infinitely large 

and this assumption is reasonable since the deflection is highly 

localized near the concentrated load location (within about 1/2 inch 

radius as seen from experiments). Recall that P = K.W so that, 

b 


K b » 16448. Ib/in 


(26) 


k b = V (E av h) = °* 029009 


(27) 
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4. Discussions of Results 


The graphite-epoxy plates and the ALS honeycomb sandwich panels were 
subjected to a dropweight impactor loading with a mass of either 2.66 lb 
or 3.9 lb. By varing the initial height of the impactor, the velocity of 
the impactor just before hitting the plate was recorded by the machine. 

The "Dynatup" IBM/PC Impact Testing System was manufactured by General 
Research Corporation, 5383 Hollister Ave., Santa Barbara, Calif. 93111 
Tel (805)-964-7724. , dated Sept. 12, 1985. 

Table 1 shows the maximum peak impact force for graphite/epoxy 
circular plates and the initial kinetic energy, assuming no energy loss. Both 

the clamped in-plane immovable boundary condition and the simply supported 

in-plane movable condition are considered and the results are tabulated 

in this table. It appears that the experimental results show that 

the plates are closer to being clamped rather than simply supported 

at the edge.. The majority of the energy was loss due to the vibration 

of the impactor system and it ranges from 70 to 75% energy loss. 

Assuming a 75% energy loss, the predicted maximinn peak impact forces 

are plotted in Figure 1 along with the experimental data. Three different 

boundary conditions are used in the experiments (i) the circular plates 

are glued to the circular blocks to avoid in-plane slipping 

(ii) no glue is applied and in-plane slipping may not be fully prevented 

(iii) the circular plate is resting on a three Inches diameter hole with 
no supporting system on top of the plate. The peak impact forces for 
the circular plates with these three boundary conditions differ from 

5 to 12% and the theoretical values agree with the experiments. 
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Figure 2 shows the peak impact force versus the initial kinetic 
energy for graphite/epoxy square plate with side being 3 inches. It can 
be seen that the measured impact forces for square plates are almost 
identical to the circular plates. This experimental observation Is 
consistent with the theoretical predictions described in section 2 
for square plates. Again, the theoretical impact forces (shown by the 
solid curves) agree with the experimental values. 


Finally, the maximum peak impact force versus the initial kinetic 
energy for honeycomb sandwich panels are shown in Figure 3. In order 
to test the validity of the theoretical model of a plate resting on 
an elastic foundation, the "loose" honeycomb sandwich panels are 
used for comparison purposes. The "loose" honeycomb sandwich panel 
consists of top and bottom face sheets and a honeycomb core but the 
face sheets are not glued to the core. Experimental data show that 
the "loose" and "bonded" honeycomb sandwich panels have approximately 
the same impact force. The theoretical impact force are tabulated in 
Table 2. Based on a 703S energy loss, the theoretical impact 

forces are plotted in solid line in Figure 3 and they agree with the 
the experimental data. The theoretical impact forces are higher than 
the experimental data since the effects of "buckling" or "crushing" of 
the core is neglected in the plate on elastic foundation model. 

In all these three figures, it can be seen that the peak impact 
force is proportional to the initial kinetic energy, at least for low. 
initial kinetic energy. 
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W/h 

P 

Clamped 

Invnovable 

(pounds) 

0/2)M I V I 2 

(ft-lb) 

P 

SS Movable 
(pounds) 

(l/2)M I Vj 2 

(ft-lb) 

0.10 

53.0350 

0.2418323 

22.8399 

0.099145 

0.20 

107.7156 

0.9534993 

46.07643 

0.392906 

0.30 

165.6876 

2.1523598 

70.01061 

0.8848849 

0.40 

228.59651 

3.876113 

95.32561 

1.5833349 

0.50 

298.08794 

6.179474 

122.1314 

2,500326 

0.60 

375.80758 

9.133731 

150.92029 

3.651588 

0.70 

463.40106 

12.826800 

182.0886 

5.056470 

0.80 

562.51402 

17.363477 

216.0332 

6.737934 

0.90 

674.79211 

22.865793 

253.1505 

8.72257 

1.00 

801.88096 

29.473412 

293.8371 

11.04062 

1.10 



338.4896 

13.72602 

1.20 



387.50479 

16.81643 

1.30 



441.27900 

20.35326 

1.40 



500.2089 

24.38172 

1.50 



564.69117 

28.95090 

1.60 



635.12231 

34.11373 

1.70 



711.89895 

39.927134 

1.80 



795.4176 

46.45196 

1.90 



886.0750 

53.753146 

2.00 



984.2677 

61.89964 


Table 1 Maximum peak impact forces and the initial kinetic Energy 
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Impact Energy (ft-lbs) 

Figure 1 Maximum Peak Impact Force versus the Initial Kinetic Energy 
for circular graphite-epoxy plates (diameter = 3 inches) 




Impact Support Comparison (Square) 
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Impact Energy (ft-lbs) 

Figure 2 Maximum Peak Impact Force versus Initial Kinetic Energy for Graphite Epoxy 
Square Plates with Side being 3 inches 




ALS Honeycomb Sandwich Panel 


W/h P (lbs) (1/2 )MjV t 2 




(ft-lbs) 

0.10 

133.2323 

0.667702 

0.20 

266.464 

2.56509 

0.30 

399.6970 

5.65577 

0.40 

532.929 

9.924744 

0.50 

666.1618 

15.3627 

0.60 

799.3941 

21.96337 

0.70 

932.626 

29.7215 

0.80 

1065.85 

38.6334 


Table 2 Predicted Peak Impact Forces using a Plate 
on a Linear Elastic Foundation Model 
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Figure 3 Maximum Peak Impact Forces versus the Initial Kinetic Energy for 
for "loose" or "bonded" honeycomb Sandwich Panels 


5. Conclusions 


The theoretical energy balance model is used to predict the maximum 
impact forces based on a given initial kinetic energy of the impactor. 

The low energy impact resistance of graphite/epoxy circular and square 
plates subjected to concentrated forces at the plate center is examined. 

A theoretical model of a plate resting on a linear elastic foundation 
is proposed for the ALS honeycomb sandwich sructure and the theoretical 
impact forces agree with the experimental data. 

Further work is being planned for the four point bending of the 
ALS sandwich panel. The sandwich panels were previously damaged and 

it is important to study the residual strength of the structure due 

to various loads including the shear loads. 
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APPENDIX A Constitutive Equations for Laminated Plates 


For a symmetrically laminated plate, there is no bending-stretching 
coupling and the membrane stress resultants(N x , Ny, N X y)are related to the 
in-plane strains of the middle surface^, ty , f ) by 


N x 


A n 

A ,2 

A 1 6 


e x 

N 

y 

| 

i 

A 12 

A 22 

A 26 

j 


N 

xy 


A 16 

A 26 

A 66 




(Al) 


Further, the bending stress resultants (M , M , M ) are related to the 

x y xy 

curvatures ( ) by 


M , 

X 

D n 

°12 


:h 


: 

M = 

y i 

D 12 

°22 

»26 

| 

ty 

(A2) 

i , 

M 

xy 

°T6 

°26 

d 66: 

) 

xy 



In the above expressions, the extensional stiffness A^ • and the bending 
stiffenss can be obtained from, 

N 

A ij = 2 Q ij (k th layer) (Z k - Z k _, ) (A3) 

k=l 

N 

3 D i j = Z Q i:j (k th layer) (Z k 3 - Z^ 3 ) (A4) 

k=l 
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Each layer is assumed to be orthotropic and the angle between the fiber 
direction and the x axis is 9. Denoting S = sin 9 and C= cos 9, the 
stress strain relations are. 



where. 


^n c4 + ( 2 ^i2 + 4( W + ^22 s4 

<T 22 = Q^S 4 + (2Q 12 + 4Q 66 ) S 2 C 2 + Q 22 C 4 

^12 ^11 + A 22 " 4 ( ^ 66 ^ ^ C 2 + A 12 ^ 4 + ^ 

{A6) 

Q 66 = ^11 + Q 22 “ 2Q 12 " 2 Q 66 } S C + % 6 (s4 + c4) 

^16 = (Q 11 ‘ Q 12 " 2( W Sc3 + (Q 12 " Q 22 + 2( W ^ 

^26 ^ Q ir Q 12 " 2C W s3q + ^U-^22 + 2 ^ Sc3 


Note that the stress strain relations in the material coordinate are, 
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The face sheet of the Advanced Launch System (ALS) structure is made 
of graphite/epoxy T300/934 material where the material parameters are, 

Ey = 20.1x10^ psi tension, Ey= 19.4 xlO E psi compression 

Ey = 1.5x1 0^ psi tension, Ey = 2.4 xlO^ psi compression 

G = 0.66 xlO E psi 

(A8) 

= 0.294 , 

density = \ = 0.057 lb-mass/in^ 


In the present analysis, the tensile values of Ey and Ey are used so that, 
E L /Ey = 13.400 

(A9) 

G L y/Ey = 0.4400 

^TL = ( V LT^ ( E j/ E L^ = 0-021940299 

From Jones (1975), eqn. 2.61, one obtains, 

(Qir ^22* Q 12^ = H/Cq) (E l » Ey, > LT Ey) 

(A10) 

Q 66 = G LT ’ c o = 1 ' V LT V TL 
Thus, in non-dimensional form, 

' q ll’ q 22 ’ q i 2 ’ q 66^ = ^ 1/E T^^ir G 22* Q 12» 

= (13.48699717, 1.006492326, 0.295908744, 0.440000) 

(All) 
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The ALS composite plate consists of 16 layers with a total thickness 


of 0.081 inch where the fiber angles are, 

(0, 45°, 90°, -45°, -45°, 90°, 45°, 0°) s 


(A12) 


Since the laminate consists of 16 equal-thickness layers, the extensional 
stiffness can be written as (h = total thickness = 0.081 inch). 



(h/4) j (0=0°) +^.(0=45°) +tf. j (0=-45°) +^(9=900} 

(A13) 


Note that, 

Q 16 (9=0) = 0, Q 16 (9=90°)= 0, 

Q 16 (9=45°) = - Q 16 (9=-45°) 


(A14) 


so that by inspection, A 16 =0 and similarly, A 2g =0. Further, q^. = ^ ./E T 


( a lT a 22 * 3 1 2 * a 66' (A 11* A 22* A 12’ A 66^D^ E T h M 

= 0/4) [q (0=0°) + 2 q. (0=45°) + ^(0=90°)] 


(A15) 


Thus, we have. 


‘11 


(1/4) [ q T1 (9=0°) + 2 q n (0=45°) + ^(9=90°)] 


a. 


1 1 


(1/4) 


L q n 


+ (2/4 ) (q^+ 2q 1? +4q 


12 H 66 


+ q 22 ) + q 22 'j = 5.729035748 
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a 22 = (l/4)[q 22 (Q=0°) + 2 q 22 (9=45°) + q 22 (0=9O 0 )} 

a 22 - (1/4) ^q 22 + (2/4) ( q n + 2q 12 +4q 66 +q 22 ) + = a„ 

a ]2 = 0/4) [q 12 (9=0°) + 2q ]2 (9=45 0 ) + q 12 (0=9O 0 ) ] 

a l2 = L^12 + ^ q ll +q 22 ~ 4q 66 +2q 12^ + q 12j = 1*813617745 

a 66 = [ q 66 + ^ 2//4 ^ q 11 +q 22" 2q 12“ 2c, 66 +2q 66^ + q 66 J = 1 * 95770900 ^ 

(A16a,b,c) 


From Jones (1975), the extensional stiffness for an isotropic homogeneous 
plate is. 


A . . = 
ij 


• E av h 


fl 


1-K 


y 

av 


av' 


, 0 

0 (1 -» jn \ 


(A17) 


Thus, the in-plane "average" Young's modulus and "average" Poisson's 


can be obtained 

from. 


5.729035784 

E I 

■ E av /('-av 2 ) 


1 .813617745 

e t 

= V av E av 

(A18) 

1 .957709001 

e t 

- E av / ?(,t 'av> 
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which implies. 


^ av = 1.813617745/5.729035784 = 0.316565966 
E flV = 7. 7323591 97x1 0 6 psi 


(A19) 


For a 16-layer laminate, the bending stiffness D^. can be computed 

using the appropriate weighting factor for each layer (Tsai and Hahn 1980, 
Table 6.6, page 234), 


3d lj * 2(1/16)3 {'69 q, j >6 = 0°) + 127?,. (0=45°) + 91? 1;j (0=90 o ) +61?^ (8.-45°) 

+ 37q. j (0=-45°)+ 19q !j (Q=90°) + 7^ .. (0=45°) + ^(9= 0°)] 

d ij = °ij /(E T h3) (A20) 

Since q 16 (0=o°)=o, q 16 (9=90°)=0, one obtains, 

3d lg = 2(1/1 6) 3 Ll34q 16 (0=45°) + 108 q 16 (0«-45 0 )] = 2(1/16) 3 26q 16 (0=45°) 

d 16 = (2/3) (1/1 6) 3 (26 /4) (q^-q^) = 0.013203659 

(A21 ) 

d 26 = d l6 = 8-01 3203659 


The remaining bending stiffness 
computed from, 


coefficients 



d 22* d 12* d 66^ can be 


d. . 

i i 



(A22) 
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so that, (d 22 /d n = 0.778487408), 


d 11 = (2/3)(1/16)3 L 170q ll + (232/4) (q 11+ 2q 12 +4q 66 + q 22 ) + HOq^ 
d n = 0.550216712 = D n /(E T h 3 ) 

d 22 = (2/3) (1/1 6) 3 [l70q 22 + (232/4) ( q 1 -j +2q-j 2 + 4q g6 + q 22 ) + 1 1 Oq -j -j j 
d 22 = 0.428336782 = D 22 /(E T h 3 ) 


d 12 = (2/3) (1/16) 3 [l 70q 12 + (232/4)(q n +q 22 - 4q 6g +2q 12 ) +110q 12 
d ]2 = 0.139277710 


(A23) 


d 66 = (2/3) ( 1 / 1 6 ) 3 [l70q 66 + (232/4) ^ ,+q^- 2q ]2 ) + HOq^'j 
d 66 = 0.151285315 


From Jones (1975), the bending stiffness for an isotropic homogeneous 


laminate is. 



Eh 3 


12(1 -v 2 ) 


1 


V 

o 


V 0 

1 0 ' 

o (1-9)72 


Therefore, by comparison. 


(424) 


0.550216712 Ej = F/ ' 1 2 ( 1 - V 2 ) j 
0.428336782 E T = F/J2(1-9 2 )] 
0. 1 3927771 0 E T = v E/ ? 12(1- 7 2 ) j 
r ‘ ' -'2853’ ? E = E/ 2<i ' 7 ) ’ 


(A25) 
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Thus, by assuming that v=0.3, one obtains, 

F = 9.0125497xlw 6 psi 
F = 7.01 61 56489x1 0 6 psi 

(A26) 

F = 7.6045629x1 0 6 psi 
F = 7.0801 5274x1 0 6 psi 

Since 77.848$, it is obvious that the bending stiffness in the 

1 -direction is not the same as that in the 2-direction. This causes 
a non-axi symmetric behavior in the deflection of a circular plate under 
a lateral concentrated force at its center. For a conservative design, 
one should use the average Young's modulus of 7.01615xl0 6 psi. 
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APPENDIX B Coefficients of Hertz Law 


According to Hertz Law, the energy due to the caontact between two 
elastic bodies is, 

E c = (2/5) P 5/3 / n 2/3 (B1 

where P is the contact force and (see Shivakumar et al . 1985b) 
n = 4 R r 1/2 / [3“ (K-j+Kg)] 

K ! = 0- Vj 2 )/( r Ej) ( B2 ) 

K 2 = (C 22 )1/2 li- (C n C 22 )1/2 + G zH 2 ~ (C 12 +G zr )2 l V2 

(2 71 )(G zr ) ^ ( C n C 22 " C 12^ 

In the above, Rj is the radius of the impactor, Ej and Vj are the Young's 
modulus and Poisson's ratio of the impactor and 


c ir C ll /E av = (E z /E av )(1_ ’ J av ^ 


C 22 C 22 /E av = ( Hi- ' 2 <i)/(l + v ) 


zr 


av 


c 10 = C 10 /E = v 

12 12' av zr 


3zr ■ V E 


av 


p * v l>- - 2 ■',} i 


$ = E /E 
av z 


(B3) 


OF PGG* 


KtvAUTf 
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Thus, equations B2 can be written in the non-dimensional form. 


n* = n/(E av h ,/2 ) = 4 (R,/h) 1/2 / [srtE^+E^)] 

V, ■ o-”i 2 >< e ™' e i >'* 


E av K 2 = 


< c 22 )’ /2 { :( c „ c 22 ) 1/2 + g 2r ] 2 - ( c 12+ 9 2r ) 2 > 1/2 


(2.t)(g 2r ) 1/2 (c 


ll c 22 


C 12 2 ) 


(B4) 


In the present analysis, we have, 

E = 7.0000 xlO 6 psi, i) , = 0.30 

ci v a v 


(B5) 


From Shivakumar et al 1985b, one can estimate the following material 
parameters , 

V = 0.060, 
zr ’ 

so that, 

i = E, w /E, = 4.117647059, 

a V Z 

p = 1 .491751492 

(c ir c 22 * C 1 2 * 9 zr^ = (°- 253597754 » 1.130491131, 0.089505090, 

0.084285714) 


G zr = 0.59xl0 6 psi , E z = 1.70xl0 6 psi 

(B6) 
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Based on the computed values of c^-j, C 22 > c-jg* 9 zr » one obtains, 


E av K 2 = le 2441 77058 


(B8) 


For the impactor, 

= 0.30, Ej = 29.0x1 0 6 psi, Rj = 0.250 in 

E K. = 0.069918413 

a v i 


(B9) 


so that, 


n* = 0.567400142 


(BIO) 


XV- 2 9 



APPENDIX C Design Guicfiine for Stacking Sequence of Quasi -Isotropic 

Laminated Plates 

From Appendix A, it was demonstrated that the bending stiffness in 
the 2-direction is only 77.848% of the bending stiffness in the 1-direction. 
The material parameters of the ALS graphite-epoxy T300/934 laminate are,. 

E l /E t = 13.400, G lt /E t = 0.44, 0.294 

and the stacking sequence is, 

(0°, 45°, 90°, -45°, -45°, 90°, 45°, 0 O ) $ 

There are sixteen layers and the total thickness is 0.081 in. From 
Tsai and Hahn. (1980 table 5.4 and table 6.6) and Jones (1975), if the 
laminate consists of "equal -thickness" lamina, the extensional stiffnesses 

(An, A^ 2 > A-j 2 » ^66’^16’ ^26^ are independent of the stacking sequence 
That is, an interchange of say any two of the above layers would not change the 
"in-plane" A^ stiffnesses and each layer carries the weighted factor 
of unity. On the other hand, the "out-of-plane" bending stiffnesses 

^11*^22* ^12* ^66’ ^16’ ^26^ depends on stacking sequence such that 
outer layers carry a larger weighted factor than the layers near the 
middle surface (halfway between the top and bottom fibers). For example 
in the small deflection bending of a beam subjected to a three point 
bending load, the top half of the beam is under compression and the 
bottom half is under tension so that the middle surface has zero in-plane 
stress and the layers near the middle surface should carry a smaller- 
weighted factor than those near the outer fibers. The weighted factors 
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for a 32 layers laminate is (Tsai and Hahn 1980), 


(721, 631, 547, 469, 397, 331, 271, 217, 169, 127, 91, 61, 37, 19,7, l) g 

in the present ALS panel which consists of 16 layers, the weighted 
factors are (there are "four" sets of the 0, 90, 45,-45 layup), 

(169, 127, 91, 61, 37, 19, 7, 1) $ 

It was demonstrated from the Utah-1 ami nate computer program that 
as one increases the number of "sets" to 16 for a 64 layer laminate, 
the bending stiffness and 0 22 are essentially the same (within 1%). 

If the number of layers is only 16, one can interchange say two layers 
so that the reatio of 0 22 to D-^ would be closer to unity. For example, 
in the above stacking sequence, one obtains, 

for the 0°, the weighted factor is 169+1 = 170 

for the 45°, the weighted factor is 127+7 = 134 

for the 90°, the weighted factor is 91+19 - 110 

for the -45? the weighted factor is 61+37 = 98 

Consider the following stacking sequence, 

(45, -45, 0, 90, 90, -45, 0, 45) $ 

one can easily show that is identical to 0 22 since 91+7 = 61+37. 
However, such stacking sequence has the drawback that bending stiffness 
in the +45 or -45 directions are not the same (169+1=170 and 127+19=146). 

A compromise stacking sequence may be. 
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(A, 6, C, D, B, 0, A, C)^ 

where the weighted factors of A, B, C, D are (169+7=176, 127+37=164, 

91+1 =92, 61+19=80) and A, B, C and D refer to 0°, 90°, 45°, -45° 
respectively or A,B,C,D, may refer to 90°, 0°, 45°, -45°, respectively, etc. 
This stacking sequence represents the "best" design since the composite 
plate behaves like an isotropic-homogeneous plate. 

Halpin (1984, section 4.6.2) suggested that the interlaminar shear 
stresses cr zx will be "significantly" lower if the +45° and -45° layers 
are separated by a 0° or 90° layer. for a quasi-isotropic layup 
involving (45, -45, 90, 0). He showed that there are only 12 distinct 
stacking possibilities. Among these 12 possibilities, six of them 
involve adjacent +45° layers and six with +45° interspered between 
0° or 90° layers. (see tables B1 and B2 from Halpin 1984). 

Thus, there is no simple solution in order to satisfy the two 
criteria (i) the laminate behaves like an isotropic-homogeneous plate 
(ii) the interlaminar stresses are minimized. It appears that 
increasing the number of "sets" to say 6 for a total of 24 layers, 
making sure that the +45° and -45° layers are separated by 0° or 90° 
layers, would be the best stacking sequence design. 
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Table Cl 


Strength 

Heirarchy 

Laminate 
Stress o v 
(KSO 

Laminate 

Interface 
Moment 
in- lb 
in 

Finite Element 
Maximum Stresses 
(KSI) 


26 5 - 

90 


-0.33 


7 

0 5 - 

0 


-1.00 

o, - -8.2 


13.5 • 

45 

i 

1 50 

t„ - 9.0 


13 5 - 

45 E 

j 

1.67 















26.5 - 

90 


-0.33 


8 

13 5 - 


J 

082 

o, - 7.4 


13.5 - 

-45 


0.98 

x„ - -9.2 


0.5 - 

0 


-0.96 



0.5 - 

0 


-0.01 


-9 

26 5 - 

90 


0.34 

Of - 7 6 


13 5 - 

45 < 

1 

085 

t„ - 9.2 


13 5 - 

-45 1 


-1.02 



0 5 - 

0~~ 


-0.01 


10 

13 5 • 

45 < 


0.15 

rt, - 10.0 


n s • 

4'. 

j 

O 115 

5. 6 1 


26 5 - 

90 C 

) 

0 98 



13 5 ^ 

45 ’ 


0 17 




^ 

) 



11 

13 5 - 

45 


0.67 

o z « 9.0 


26 5 - 

90 


1.02 

- -7.7 


0.5 - 

o i 

1 

1 02 



13 5 - 


■N 

0.17 


1? 

13 5 • 

45 ( 


0 67 

o, = 10.9 


05 • 

O 


1 34 

5. 1 ? 


26 5 - 

90 r 


1.67 




G 

] 




Q mj ■ i morn 'r 
Q maximum o t 

T/300 5208 Graphite Epoxy 
i, - 0 5% 


Normal Stress and Interlaminar Stress for adjacent +45° 
Quasi-isotropic laminate 
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T ( 300^5206 Graphite- Epoxy 
i a '0 5% 



Stress and Interlaminar Stress 
with either a 0° or 90° Quasi - 


for Interspersed 
isotropic Laminate 


Table C2 


Appendix D Equilibrium and Compaitbil ity equation for Circular Plates 

and Rectangular Plates 

From Chia(l 980, section 3.2), the nonlinear equilibrium and compatibility 

equations (written in terms of the out-of-plane displacement W and a stress 
function F) are, respectively, assuming axi-symnetric deflection, 

(D)(R)[(1/R)(RW, r ), r ] , r 

r C(i/r)(rf, r ), r 3, r 


s=R 

i 

s=0 


sq ds + F, D W, 


(Dl) 


R *R 


(-Eh/2) (W, R )‘ 


(D2) 


where 


N r = (1 /R)F, r 


N = F 
'9 r, RR 


(03) 


N Q = 0 
r9 


In the above expression, R is the radial coordinate, E is Young's modulus, 
h is the thickness, q is the applied stress, and N r , Ng are the membrane 
stress resultants. Assuming that the concentrated load can be represented 
by an applied stress over the contact area of radius a con t ac t’ anc * dropping 
the nonlinear term F, R and W, R in the equilibrium equation for small deflection, 
one can show that the exact deflection is of the form. 


W = W-j (1 + Cl r 2 + c 2 r 2+<1 ) 


( 04 ) 


where : is arbitrarily close to zero 
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The governing nonlinear equilibrium and compatibility equations "for 
symmetrically laminated rectangular plates are, respectively, (Hui 1985a, b) 


l d *(w) F »YY W ’XX + F 'XX W *YY " 2F, XY W, XY 


(D5) 


L A*( p ) = 


(W, xy ) - W, xx W, Yy 


(D6) 


where W is the out-of-plane deflection, F is the stress function,.'. 

X and Y are the in-plane coordinates and L Q *( ) and l ft *( ) are the 
differential operators defined by, 

l a *( ) = A 22 ( ), mx + (2A 12 +A 66 ) ( ), XXYY + A-j i ( ), yYYY 
l d *( ) = D n( )» xxxx + ( 2D i2 +4D 66^ )*XXYY + D 22^'YYYY 

★ 

The A.j and are the material parameters defined by Jones (1975). 
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Finally, an extension of the work by Eringen(1953) to laminated plates 
was presented by Sun and Chattopadhyay (1975) based on the summation 
of the various deflection modes corresponding to different frequencies 
of an anisotropic rectangular palte. This method was used by Chou 
and Mortimer (1976) who presented a computer code to predict the contact 
force, deflection and bending strains of an anisotropic plate. Using 
this code, the predicated strains were in good agreement with the experimental 
data as reported by Dobyns and Porter (1981). 

The present work deals with the prediction of the peak impact 
force due to low energy impact on graphite/epoxy circular a<id square plates 
and ALS honeycomb sandwich panels. Using the energy balance method, 
the initial kinetic energy is dissipated in terms of the bending energy 
of the plate, the membrane energy(due to the stretching of the mid-surface 
of the plate in finite deflection), the Hertzian contact energy (due to the 
imbeddment of the impact in the plate) and the energy loss (due to the 
vibration of the impactor punching system). The "average" Young's modulus 
and "average" Poisson's ratio as well as the estimated out-of-plane 
material parameters (see Shivakumar and Crews 1982 andKriz and Stinchcomb 
1979) are reported based on the ALS 16-layer quasi-isotropic layup. 

The Hertzian parameters are then computed in Appendix B. The theoretical 
contact forces are in good agreement with the experiments. Some design 
guidlines for various stacking sequence are presented in Appendix C. 
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